Screening effects on ^5*0 pairing in neutron matter 
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The ^So superfluidity of neutron matter is studied in the framework of the generahzed Gorkov 
equation. The vertex corrections to the pairing interaction and the self-energy corrections are 
introduced and approximated on the same footing in the gap equation. A suppression of the pairing 
gap by more than 50% with respect to the BCS prediction is found, which deeply changes the 
scenario for the dynamical and thermal evolution of neutron stars. 
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Neutron superfluidity in neutron star matter though an 
old subject, is still of great actuality and vividly debated 
(for a review, see The reason for this stems from the 
fact that superfluidity is an extraordinarily subtle process 
when it comes to quantitative predictions starting from 
the bare NN interaction. On the other hand for neutron 
stars, such quantitative predictions are necessary, since 
the manifestation of the superfluidity are only rather in- 
direct through glitches and relaxation phenomena and 
cooling rates. One, therefore, lacks direct experimental 
information on the magnitude of neutron pairing. On the 
other hand, there is no doubt that the dynamics and the 
thermodynamics of neutron stars are strongly influenced 
by the superfluid character of neutron matter and it is 
therefore important to get pairing properties of neutron- 
star matter under better control from the microscopic 
point of view. Of course, neutron matter superfluidity 
is not completely decoupled from the one prevailing in 
finite nuclei. Though we have experimental information 
for these objects, also in this case the fully microscopic 
explanation of the observed phenomena is far from be- 
ing completely settled. It may be argued that in finite 
microscopic systems, where the surface plays a very im- 
portant role, the situation can be quite different from 
the homogenous case. However, like for other quantities 
of nuclear physics it should be possible to disentangle 
volume and surface effects also for pairing properties in 
finite nuclei and it is therefore our belief that the topic 
of superfluidity in neutron matter, nuclear matter and 
finite nuclei should be studied in an interrelated way. 

In this work we again concentrate on neutron mat- 
ter in pursuing previous studies. However, it is planned 
to parallel this work for nuclear matter in the near fu- 
ture. In the past we have mainly been concentrating on 
the infiuence of either dynamic self-energy corrections 
(see also 0) or vertex corrections to the neutron matter 
pairing problem All investigations in this direction 
invariably led to the conclusion that dynamic self-energy 
corrections yield a quite strong reduction of pairing 
in neutron matter. However, to be consistent, self-energy 
corrections have to be followed by vertex corrections on 



the same footing. This is the objective of the present 
investigation. It will be seen that the vertex corrections 
have tendency to further reduce the gap but to a lesser 
extent than that is the case from self-energy corrections. 
Our approach is based on the Gorkov Green's function 
formalism where we develop systematically self-energy 
and effective pairing interaction to lowest order in the 
particle-hole bubble insertion. Though intuitively quite 
reasonable and physically motivated, we should mention 
that this is neither based on an expansion in a small 
parameter nor does it rely on some variational princi- 
ple. None the less our results will be quite comparable 
to those of other works, notably to those of Clark et al. 
Q who based their investigation on the correlated basis 
function approach which, to a large extent, is variational. 

The superfiuid phase of a homogenous system of 
fermions is described by the pairing field Afe(cj), which 
is the solution of the generalized gap equation 
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Here V is the sum of all irreducible nucleon-nucleon (NN) 
interaction terms and Fk{Lu) the anomalous propagator 
JJ- In most pairing calculations ^ the effects of the 
medium polarization has been included in the self-energy 
0, Qj E^'^d not in the pairing potential. A more gen- 
eral study requires the medium corrections to be treated 
on equal footing also in the vertex corrections as well 
as in the self-energy. This is the main concern of the 
present letter. Accordingly, the expansion of the interac- 
tion block V and the self-energy S have both been trun- 
cated to second order of the interaction. The correspond- 
ing diagrams are depicted in Fig. 1. 

The limitation to lowest order bubble insertion may 
seem as a strong restriction because resummation of the 
bubble series into RPA can have an important infiuence. 
However, in the exploratory work where self-energy and 
vertex corrections are treated consistently for the first 
time, we think that higher order effects unnecessarily 
complicate the approach and that the lowest order effects 
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in the density should at least give the correct tendency 
even up to densities around saturation. 

The matrix elements of the bare interaction V exhibit 
hard-core divergences which have to be removed by dress- 
ing them with the short-range particle-particle correla- 
tions (ladder diagrams). This can be done either in a 
microscopic approach by replacing V with the G-matrix 
or in a semi-phenomenological approach by replacing V 
with an effective interaction. 

In the calculations of all diagrams, the superfluid prop- 
agators have been replaced by the normal phase propa- 
gators. This amounts to neglecting second-order correc- 
tions in the gap, which are actually negligible as shown 
in a fully self-consistent calculation discussed below. 
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FIG. 1: The diagrams of NN interaction and self-energy dis- 
cussed in the text. The exchange terms are understood. 

In general, the interaction as well as the self-energy are 
complex quantities which implies the energy gap to be 
a complex function. The introduction of the imaginary 
part of the self-energy amounts to taking into account the 
effects of the quasi-particle spectral function, which have 
been studied elsewhere 10, JTj. The complex nature of 
the potential is due to finite time propagation and decay 
of processes shown in Fig. 1. In the present investigation 
V will be assumed to be a real function what implies the 
gap function to be real also. 
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FIG. 2: (a) Rearrangement contributions to the self-energy, 
where k is fixed to kp- (b) The HF mean field is plotted vs 
momentum k at kp = 0.5, 0.8, 1.1 fm~^. 

For the present calculation we adopted the semi- 
phenomenological approach and chose the Gogny force 



Dl as effective interaction at each coupling vertex shown 
in Fig. 1. The replacement of the vertices in Fig. 1 by the 
Gogny force which can be considered as a phenomenolog- 
ical representation of a G-matrix may seem unjustified 
for the lowest order term in Fig. 1(a), since perturbation 
theory tells us that it should be the bare NN interaction. 
However, it is well known that limiting the fc-space, 
the bare interaction has to been replaced by an effective 
one and choosing the cut off to be at kp the effective 
force in Fig. 1(a) can be shown to be again equivalent to 
the G-matrix. In neutron matter only the density inde- 
pendent part of the Gogny force survives with a range in 
/c-space of about kp at the saturation density. Therefore, 
from this point of view, it may not be unreasonable to 
take the Gogny force even for the lowest order term in 
Fig. 1(a), which is further backed by the fact that the 
gap with the Gogny force (without polarization terms) is 
quite close to the one calculated from a bare NN force. 
In any case the qualitative influence of the polarization 
on the lowest order solution of the gap equation, given 
by retaining only graphs Fig. 1(a) and Fig. 1(c), should 
not depend very much whether we use in Fig. 1(a) the 
bare interaction or the Gogny force. 

The self-energy T,k(uj) in neutron matter has been cal- 
culated within the above described approximation. The 
first-order contribution is reported in Fig. 2(b) for three 
values of the nuclear matter density in the range where 
the pairing is expected to be largest. The second-order 
contributions T,pp and Sj^^^ (graph b) and c)) in Fig. 1, 
respectively) are depicted in Fig. 2(a) as a function of 
the energy. Only the energy dependence of of these terms 
will be discussed, since the Gogny force implicitly con- 
tains the static part, which will be removed from the gap 
equation as described later on. 

(2) . . 

SJj/j exhibits a pronounced maximum in the vicinity of 
the Fermi energy due to the high probability amplitude 
for particle- hole excitations near ep- It is in very good 
agreement with the results obtained from Brueckner- 
Hartree-Fock (BHF) calculation with G-matrix 2]. 

The second-order potential is given by the one-bubble 
exchange term (plotted in Fig. 1(b), which is the first 
one of the ring diagram series. Physically it represents 
the screening to the pairing due to the medium polariza- 
tion. Again we take the Gogny force for all vertices in 
Fig. 1 (a,b) . Our prediction for Vj.^,, {ui) at three typical 
densities is reported in Fig. 3. We plot the symmetric 
part 
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;[V(^,c^^) + V(^,-4)], (2) 



which is the only one relevant for the pairing gap. The 
strength of V is concentrated around the Fermi energy 
(lu = 0) with a peak value at k = k' — kp and a width in- 
creasing with the density. Its uj dependence is shaped by 
the polarization part, i.e., Lindhard functions [l2j | which 
at cij = (static limit) is repulsive at any momentum 
and density, but it becomes attractive for |a;| ^ ep. One 
therefore expects a reduction of the gap due to screening. 
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FIG. 3: Screening potential vs energy at kp = 0.5, 0.8, 1.1 
fm~^, separately. 
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The two poles are located close to the real axis on oppo- 
site sides of the imaginary axis. Leaving aside a general 
integration of the gap equation, we adopt the pole ap- 
proximation relying on replacing the full propagator by 
its pole part: 



ZfcAfe(w) 



1 



1 



In general the residue Zk at the poles is defined as 



Zk=\l- 
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calculated at the Fermi surface. In the calculations we 
took the limit of for A ^ 0, which corresponds to the 
quasi-particle strength. The factor Z is keeping the full 
dynamical dependence of the self-energy reported in Fig. 
2. Afterwards we may subtract the static part E(^)(w) 
from the self-energy. In the calculation of the gap the 
inclusion of the latter only brings a variation of less than 
one percent since the gap is not sensitive to the static self- 
energy far from ep- Inserting Eq. (8) for the anomalous 
propagator into the gap equation, after u integration we 
obtain 



k' dk'Vk,k'{^ - ^k')- 



Zk'Ak'{oJk') 



(10) 

Notice that the reason why only the w-even part of 
the interaction contributes to the integral can be traced 
to time reversal invariance of the superfluid ground state 
for which the anomalous propagator as well as the gap 
function are even functions of uj. 



The gap equation, Eq. (1), is to be coupled to the 
closure equation for the Green's function 



(3) 



fixing the chemical potential in the superfluid phase. 
The anomalous propagator 



Fk{uj) 



[huj - ek(uj)] [hoj + efc(-w)] ' 



(4) 



has two poles, symmetric with respect to the imaginary 
axis of the complex w-plane 3], which are the roots of 
the equation 

±ujk ^ £ki±(^k)- (5) 
The quasi-particle energy is given by 



Skicu) = J:zico) + W el + + AUco), (6) 
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FIG. 4: Energy gap in the present approximation. For com- 
parison the prediction from the pure BCS model (dotted line) 
and from BCS plus self-energy effects (dashed line) are plot- 
ted. 

The remarkable advantage of this approximation is 
that the gap depends only parametrically on lo and its 
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energy dependence is only related to the energy depen- 
dence of the interaction. The on-shell gap Afe(ajfc) fulfills 
the equation 

(11) 

equivalent to the gap equation in the static limit. 

The approximate version of the gap equation, Eq. (10), 
has been solved using as input the self-energy and pair- 
ing potential discussed in the previous section. The fo- 
cus was in the ^Sq neutron- neutron pairing for neutron 
matter, which is by far the most important component 
of pairing as to possible implications in nuclear systems. 
The energy gap at the Fermi surface (fc = hp and w = 0) 
as a function of fc^ is reported in Fig. 4. The domain 
of existence of the superfluid state is mainly at low den- 
sities with a peak value of about 1.4 MeV at hp = 0.6 
fm^^. In the same figure we also report the result in 
the BCS limit [for a review see U] (neither self-energy 
effects nor screening) and the result with self-energy ef- 
fects but without screening 0. From the comparison of 
the three predictions one sees that the main suppression 
of A is due to the strong g.s. correlations which lead 
to a Z- factor much less than unit. However screening of 
the pairing interaction produces an additional suppres- 
sion. It has to be noticed that the screening potential 
also shifts the peak value of the gap to lower density, 
where the suppression is less sizeable. This implies that 
the medium effects are not expected to reduce dramat- 
ically the pairing at the nuclear surface and the feature 
of pairing as surface effect is still confirmed by our fully 
consistent calculation of pairing in neutron matter. As 
final result, obtained solving Eq. (10), we report in Fig. 
5 the gap as a function of the energy for kp = 0.6 fm^^. 
We address here the relevance of such a feature for the 
study of pair correlations in dynamical processes such as 
the expanding and disassembling phases of the fragmen- 
tation events in heavy-ion collisions jl7j |. 
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FIG. 5: Energy gap as a function of u at kp = 0.6 fm~^. 

Screening effects on the pairing interaction have also 
been studied in different contexts. One of the earli- 



est calculations has been performed in the framework of 
the second-order correlated basis perturbation theory , 
where a pairing suppression by a factor 4 is predicted. 
In that work the momentum dependence is mostly ne- 
glected that amounts to overestimate the suppression. 
Closer to the present approach is the polarization poten- 
tial calculated from the induced interaction theory 
[1^ , which gives substantially the same gap when treated 
within the Landau parameter approximation fl^ . Fi- 
nally, we should mention a parallel study on finite nuclei 
|l6l | , where the polarization potential is given by the cou- 
pling to surface vibrations. While the self-energy plays 
the same role as in neutron matter, the phonon exchange 
produces an enhancement of the pairing at variance with 
the dominant repulsive effect of the spin fluctuations in 
neutron matter. 

In conclusion, this work constitutes a continuation of a 
previous one(Ref. '^) where we investigated self-energy 
effects on the pairing gap in infinite neutron matter. Here 
we treated consistently the additional inclusion of ver- 
tex corrections. On the same footing we considered the 
lowest order particle hole polarization bubble both in 
the self-energy and in the screening of the pairing force. 
Instead of the G-matrix we used the phenomenological 
Gogny force at all coupling vertices. We verified that this 
replacement has only very little influence on the numer- 
ical results. The screened pairing interaction is in prin- 
ciple energy dependent but in the quasi-particle approx- 
imation used here, this dependence on energy becomes 
only a parametrical one which greatly facilitates the nu- 
merical task of solving the gap equation. The outcome 
of the inclusion of vertex corrections is that the gap as 
a function of kp maintains approximately its bell-shaped 
form, but with respect to the self-energy corrections only, 
a further substantial reduction of the gap is induced by 
screening the bare NN interaction in the gap equation. 
Therefore with respect to the lowest order approach, i.e., 
without any polarization effects, that is with bare inter- 
action and k-mass only, the gap at its maximum is now 
reduced by about 50%. This strong reduction is a com- 
mon feature of all previous calculations and in this sense 
our investigation is a confirmation of what has been found 
by other authors earlier, even though the approaches dif- 
fer in detail. Nonetheless this strong reduction of the 
pairing due to the polarization remains intriguing. If the 
same situation should prevail in nuclear matter an esti- 
mate via the local density approximation |^ would lead 
to a by far too low value of the gap in finite nuclei. How- 
ever the influence of polarization terms in nuclear matter, 
due to the different quantum numbers involved, may be 
quite different from the one in neutron matter. It may be 
an extremely interesting problem for studies in the near 
future to see whether the polarization corrections in the 
different channels, i.e., n-n pairing in neutron matter, n-n 
pairing in nuclear matter and n-p pairing in nuclear mat- 
ter, can, at least qualitatively, explain the experimental 
finding known from finite nuclei. 
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